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By Wertheim-method the exact solution of the Percus-Yevick integral equation for a sys- 
tem of particles with the "repulsive step potential" interacting ("collapsing" hard spheres) is 
obtained. On the basis of this solution the state equation for the "repulsive step potential" is 
built and determined, that the Percus-Yevick equation does not show phase transition for "col- 
lapsing" hard spheres. 
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PACS: 61.20.Ja,64.70.Ja 

In 1963 Wertheim and Thiele indepen- 
dently gained the exact analytical solution of 
the Percus-Yevick integral equation for hard 
spheres P-|S1- This solution is so far the only 
strict analytical one of non-linear integral 
equation for pair distribution function. The 
given work shows the possibility of Wertheim 
method to express the solution of the Percus- 
Yevick equation in closed analytical form for 
more complicated particles interaction poten- 
tial — repulsive step potential j4j 



V{r) = 




(1) 



where Vq — positive constant, r — particles 
range ( "collapsing" hard spheres) . 

Such a potential used to have a great 
application in modeling phase transition in 
liquid under the high pressure, isostructural 
phase transitions in crystals, transformations 
in colloid systems, etc. by means of molecular 
dynamics within the framework of thermody- 
namic perturbation theory [U |3] . 
Let's consider the Percus-Yevick equation 



n2{r) = 1 — n 



[n2{f — s) — 1] ds 



where n2{r) — the pair distribution function, 

/?=-—, n — the particle density. Moving to 

kT 
bipolar coordinates and integrating for angle 

variable for the "repulsive step potential" we 

gain that 



r+s 

h{r) = Ar- 2-Kn I h{s)ds f h{t)e~^^^''^dt 
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(3) 
where 

h{r) = rn2{r)e^^^'^' = 

-rC{r), r < a 

-rC{r) (4) 

, a < r < b 



1 - e-/5^o 

712 (r), r > b 

C{r) — direct correlation function. In the 
approach of Percus-Yevick 

C{r) = (1 - e^^("))n2(r) 

n2{r) = at r < a, C{r) = at r > 6, and 
e(t-a)e(6-t) + 0(t-6); 
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Q{x) — Haviside step function; 

a 

A = l + 47m sh{s)ds+ 



b 

+47rn(l -e"^^°) sh{s)ds 



where 

Co = 2^ne~^^>(AiA;2 + koij6i - k)) 

irnXi 



a 



12 



We take the one-side Laplace transform for 

oo 

(3) L{h{r)) = J h{r)e~^^dr and change the 



order of integration for r and t, we finally ob- 
tain: 



A + 'yz5{z) 



-Liz) 



and the constants Xi, X2,ko,ki,k2, Iq and 5i as 
density, temperature and potential parame- 
tres Vo,a, 6 functions can be obtained from 
the system of equations 

a 

Xi= A = l + Aim I rh{r)dr - 47to(1 - e~'^^'')ki 
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2lTn r . X / X 

1 [Liz) - Li-z) 



(5) 



As = 27(5i - 2/o 
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where 

ijiz) = L(rn2(r)) = Giz) + e~^^°Kiz) 

Liz) = Li-rCir)) = F(^) + (1 - e-'^^'')Kiz) 

a 


b 

Kiz) = [ his)e~''ds 



X 
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I r^hir)dr + (1 - e"^^') / r^hir)dr 
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r"'hir)dr, m = 0,1,2, 



Iq= I hir)dr + (1 - e~^^o)A;o 
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oo 

Giz) = f his)e-''ds 

b 

6iz) = aiz) — ai—z) 
7 = 27rne-^^«(l-e-'^^°) 

For further investigation we lead in the fol- 
lowing function 



Si 




his)hit)it - s)dtds 



(8) 



Hiz) = z^i/jiz) 



A + 'jzSiz) 



-Li-z] 



(6) 



Discussions like in |2j show that 

Hiz) = Ai + X2z\ 

where Ai , A2 — are contstants. Not includ- 
ing ipiz)Li—z) from (5) and (6) and turning 
the Laplace transform into the r <h area, we 
gain the explicit expression for /i(r) 

hir) = -(Co + Cir + C2r2 + C3/) (7) 



Inserting (7) into (5), we gain the Laplace im- 
age form for m2(r). 

The system of "collapsing" hard spheres 
state equation can be shown as 

— 7^ = 1 777; rn2ir) ( — — | dr = 

nkT 6kT J ^ ' \dr ) 

= 1 + ^ [e-Py^a'ria) + (1 - e-^^°)6V(6)] 

(9) 

where rir) = /i(r)r^^ and inverse isothermic 
compressibility 



\d^)^kT 

= \-n fcir)dr = Xiin,T) 



(10) 



It can be seen from (10) and (8) that if H > References 



dP 

dn 



>0 
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i.e. the Van der Waals loop is absent in the 
isotherm. This resuh coincides with those 
ones wich were taken in numerical analysis 
of the state equation (9). 
Thus, the Percus-Yevick equation for the sys- 
tem of "collapsing" hard spheres allow the so- 
lution in closed analytical form going into the 
Wertheim-Thiele classical solution for hard 
spheres when a = b. As in the case of hard 
spheres the Percus-Yevick solution doesnt 
show the phase transition in the system of 
"collapsing" hard spheres. 
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